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i Three-dimensional icosahedral random tilings with rhombohedral cells are stud- 

ied in the semi-entropic model. We introduce a global energy measure defined by 
the variance of the quasilattice points in the orthogonal space. The internal en- 
ergy, the specific heat, the configuration entropy and the sheet magnetization (as 
defined by Dotera and Steinhardt [Phys. Rev. Lett. 72 (1994) 1670]) have been 
calculated. The specific heat shows a Schottky anomaly which might indicate a 
phase transition from an ordered quasicrystal to a random tiling. But the diver- 
gence with the sample size as well as the divergence of the susceptibility are too 
small to pinpoint the phase transition conclusively. The self-diffusion coefficients 
closely follow an Arrhenius law, but show plateaus at intermediate temperature 
£ — ' ranges which are explained by energy barriers between different tiling configura- 

tions due to the harmonic energy measure. There exists a correlation between the 
temperature behavior of the self-diffusion coefficient and the frequency of vertices 
which are able to flip (simpletons). Furthermore we demonstrate that the radial 
distribution function and the radial structure factor only depend slightly on the 
random tiling configuration. Hence, radially symmetric pair potentials lead to an 
energetical equidistribution of all configurations of a canonical random tiling en- 
semble and do not enforce matching rules. 
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I. INTRODUCTION 



C . 

The stability of quasicrystals has been a riddle since their discovery in 1982. Do they form 
stable or metastable states? Are they stable only at high temperatures? Is the stability due to 
energetic or due to entropic reasons? The random tiling model of quasicrystals is an abstraction 
which deals with rigid tiles, thereby neglecting thermal fluctuations and the phonon degrees of 
freedom. The only dynamic process is the local rearrangement of tiles, called "flips" , "umklapps" , 
"zipper" moves, depending on the type of tiles and symmetries. The random tiling model has 
been proposed by Elser Q and has been studied intensively in recent years. The first author 
dealing with random rhombohedral tilings in three dimensions was Tang [|| who was interested 
in diffuse scattering and phason elastic constants. Strandburg j^] calculated the configurational 
entropy. The state of the art of random tilings was reviewed by Henley || . Ebinger || studied 
random tilings at infinite temperatures. Dotera and Steinhardt || introduced the concept of sheet 
magnetization as an order parameter to describe the randomness. Interest in random tilings was 
renewed by Kalugin and Katz [nlthrough the new process of flip diffusion. This property has been 
studied by Jaric and S0rensen||[| at infinite temperature, and by Joseph [jl0| and Gahler jll[] 
at finite temperatures. Meanwhile for the energetic interaction of the tiles many different sets of 
matching rules exist. 

In the present work we first deal with quantities that allow to characterize a random tiling. We 
are especially interested in the atomic plane and line structures important for example for ion 



channeling (Sec. III). Energy measures (Sec. IV) are discussed and a new method for generating 



non-equilibrium quasicrystals (Sec. |V|) is established. The possibility of a random tiling equilibrium 



phase transition is studied in Sec. KTf. Structure functions are introduced in Sec. VIA, and results 



concerning a special "harmonic energy measure" as defined in Sec. IV are presented in Sec. VIII 



Conclusions are drawn in Sec. IX 
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II. DEFINITIONS 



Quasicrystals can be described as cuts through higher-dimensional periodic crystals. The ad- 
ditional dimensions are addressed as orthogonal space K 1 -. Quasilattice points x in the physical 
space can uniquely be lifted to the higher-dimensional space X = II|| x x and then be projected 
onto the orthogonal space by y = n_j_(X) (II|| + = id in a proper normalization). The in- 
finitely extended quasilattice thus is contracted into a finite volume called "acceptance domain" or 
"atomic hypersurface" . The whole procedure of lifting and projecting into the orthogonal space is 
called "dualization" y = II j_ o 1177 x of the quasilattice. The higher-dimensional embedding gen- 
erates new degrees of freedom in addition to the ordinary phonons in periodic crystals, denoted 
"phasons" . 

The icosahedral quasilattice lifted into higher-dimensional space forms a three-dimensional hy- 
persurface, called "de-Bruijn- hypersurface" or "Weiringia roof" The hypersurface fluctuates 
around an average hyperplane h(x) = ho + £ g iobai~x- The constant quantity 

£ g iobai '■= V|| ® h(x) = const 

is called "global phason strain" and describes the deviation of the slope of the average hyperplane 
h(x) from the slope of the physical space. In the case of exact icosahedral symmetry the average 
hyperplane is running parallel to the physical space and thus e g i bai = 0. 

Long-wavelength deviations from the hypersurface are denoted as "phason fluctuations" and are 
described in a continuum picture by a phason strain tensor 

e(x) := V|| ® h(x) . 

According to Henley the fluctuations are governed by a free energy depending quadratically 
on the phason strain: 

F = J d d * trace (e(x) <g) e(x) T ) . 

The dimension d of physical space is equal to three for icosahedral quasicrystals. 

The present work deals with the three-dimensional Ammann-Kramer- Penrose tiling and its ran- 
domizations. This tiling consists of two different elementary cells, the "oblate" and "prolate" 
rhombohedron. Both of them appear in ten different orientations in the tiling. The six ratio- 
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nal linearly independent unit edge vectors tj may be defined as t a — -^-(cos §7ra, sin §7ra, \) 
(ct=0,. . .,4) and t5 = (0,0,l). The atomic hypersurface is a rhombic triacontahedron of icosahedral 
symmetry. The tiling exhibits, among others, a vertex representing lattice points where two pro- 
late and two oblate rhombohedra meet. It is denoted "simpleton vertex" . The outer shape of the 
cells, which touch this vertex, is a rhombic dodecahedron. If the dodecahedron is kept fixed, two 
possibilities exist to fill it with tiles. The exchange of one configuration by the other is a "flip" • 
The lattice point jumps a distance which is 0.650 of the edge length of the rhombohedra (Fig. I]). 
Any non-ordered space filling arrangement of rhombohedra without gaps or overlaps is called a 
"random tiling" . 

III. RANDOM TILING CHARACTERIZATION 
A. Variance 

Any random tiling may be characterized by the mean square deviation of the point distribution 
from the center of mass in the orthogonal space. The variance is defined by 



2 



^ N N * ^ N 

nJ2 Vi~ = iV^> J 



yo 



(i) 



N is the number of quasilattice sites and yo is the average of the position vectors yj of all dual 
quasilattice sites. The bigger the phason fluctuations are the larger is the disorder and therefore 
the random tiling parameter f2. The averaged value in the pure entropic random tiling ensemble 
amounts to 1.73±0.01 @. 



B. Vertex frequency 

The lattice points of a quasicrystal can be characterized by their local environment. The number 
of rhombohedra adjacent to a lattice point in an icosahedral quasilattice varies from 4 to 20, 
the number of edges from 4 to 12. There are 24 "canonical" or "allowed" vertices, but 5450 
vertices may occur in a random tiling. Vertices where two identical rhombohedra meet are called 
"crystallographic" . Without them the number of vertices is reduced to 360 |L4j . 

In the two-dimensional Penrose tiling each flip generates non-canonical vertices. There the 
knowledge of the vertices suffices to distinguish between ideal and random tilings. This statement 
is not valid in three dimensions. There are some flips possible in the rhombohedra tiling which only 
change the frequency of vertices without introducing forbidden ones. If the degree of randomization 
has reached a certain level the number of forbidden vertices is going to rise rapidly. 



C. Sheet magnetization 

The frequency of the simpleton is 23.61% in the ideal rhombohedra tiling. The simpletons are 
arranged in two-dimensional layers perpendicular to the two- fold symmetry axis. Due to the two 
possibilities to pack the rhombohedra into the simpleton there are two positions for the internal 
lattice point. The positions may be called "up" and "down" and attributed a spin value S = ±1. 
In the ideal rhombohedra tiling all the spins in a certain layer carry spin +1 or -1. Dotera and 
Steinhardt [|| have defined a sheet magnetization of value 1 using a proper summation rule. In a 
random tiling the sheets also exist, but the magnetization is reduced since the spins in one layer 
are not all aligned. The collection of all sheets can therefore be regarded as a two-dimensional 
Ising model with a variable number of spins. There is, however, no interaction leading to an Ising 
Hamiltonian. The sheet magnetization simply is being used as an order parameter which may 
indicate an order-disorder transition. The number of spins is not constant as in a standard Ising 
system. The susceptibility is given by 

M^\_/M_\ 2 \ 
K/ \Nn) ) ■ 

Nd is the number of the dodecahedra and therefore the number of the "spins" in the patch. This 
equation replaces the more cumbersome definition used in H. 




D. Tube and slab structures 



In the three-dimensional icosahedral rhombohedra tiling special tube structures, denoted 
"worms", exist along two- fold axes (Fig. ||). They are formed of stacks of rhombohedra which 
have one type of face in common, perpendicular to the two-fold axis. In addition to the tube 
structures there are slab structures perpendicular to a five-fold direction. The slabs are formed 
by rhombohedra which have one type of edge in common (parallel to the five- fold axis). The pro- 
jection of the slab structures onto the plane perpendicular to the appropriate five-fold direction 
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yields generalized Penrose tilings in two dimensions. These tube and slab structures exist in ideal 
tilings as well as in random tilings. 



E. Plane and chain density diagrams 



Atoms in crystals are arranged in lattice planes characterized by a normal vector of rational 
components (the Miller indices), and in lattice chains with a rational direction vector. This fact 
also holds for quasicrystals. The lattice planes or chains of a fixed direction in a crystal all are 
equivalent to one another in the simple cases (or they fall into a small number of classes). This 
is not true for a quasicrystal: For example, the density of atoms in a certain plane does not only 
depend on the orientation of the plane but also on its position. The basic theory can be found in 



Ref. [p.5[ . The plane and chain structures have been studied in Ref. [p.6 17\. 

As described in Sec. □] it is possible to lift the quasilattice points into the higher-dimensional 
space and to project them onto the orthogonal space. Points in the same physical plane are 
situated in the same plane in the orthogonal space p5| , and points on one chain in the physical 
space will lie on one chain in the orthogonal space, in the ideal tiling within the acceptance domain, 
respectively. 

To describe densities of atoms in the planes of our (primitively decorated) quasicrystal, we use 
the function A s (£^). The index s shall indicate a two-, three- or five-fold symmetry direction. 
is the coordinate of position of a plane in the direction parallel to its normal and is the dual of 
the height £" of the selected plane in the physical space (with respect to the origin for example). 
yl s (^) is proportional to the area of the intersection plane at £ through the triacontahcdral 
acceptance domain in the ideal tiling, since the density of points is 1 in the domain and outside. 

For random tilings we have to use another definition. The area of A s (£-*-) at the position £ x in 
E 1 - is proportional to the plane decoration density Pf,s(^[S, ± }) ='■ Pf,s(^) of the corresponding 
plane = const in physical space. The index F denotes plane densities. 

In the case of an ideal quasicrystal, the values of A s (£ ) are between and a maximal value 
for five- fold and three- fold directions (Fig. ||). In the five- fold case the function has a constant 
plateau since the corners of the intersection polygon are formed by edges of the triacontahedron 
perpendicular to the plane. In the two- fold case the minimal value of A s (^ x ) is greater than zero, 
since the smallest intersection polygons coincide with the surface-rhombi of the triacontahedron 
parallel to the intersection plane. To obtain the average density of atoms in a plane one has to 
calculate the mean value of the area of the intersection polygons. It is equivalent to the volume 
of the triacontahedron divided by the diameter A s of the triacontahedron along the direction s. 
Then the mean value of the plane density is: 

t Vxria 
PF,s OC A s = — . 



For channeling experiments one is interested in the frequency H(pp s ) of planes perpendicular 
to s with a plane density pp, s m the interval 8pp s . Wherever the intersection function A s (£ ) oc 
pF <s {i\\) is flat, there are many £" [£"""] values for pp, s inside of 5pF, s . If the slope is large, the 
number of planes with this density is small. The densities H(pp :S ) are obtained from the values 
Pf,s by a subdivision of the pf,s axis into intervals 8pF,s of equal length and calculating for each 
interval 8pF,s the length of the inverse image on the ^-axis. Since the triacontahedron is a 
convex polyhedron and centrosymmetric, the function t4 s (^) is monotonously growing in one 
half of the interval A s and monotonously decreasing in the other. Therefore one can calculate the 
first derivative of the inverse image (that means SpF.s — * 0) in one half and set it proportional to 
H(pf. s ) which hence is ^'(Ag). Since ^4 S (C X ) has horizontal lines and kinks, H(pp s ) possesses 
poles and jumps. 

Similarly as the plane structures one also finds chains in the quasicrystal. The direction again 
is denoted by s, the position in the perpendicular space now is characterized by a two-dimensional 
parameter {^"SC^} = const with its counterpart in the physical space at {77", C"} — const. The 
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chain density pi tS {v ■> C ) a t [v i C ] is proportional to the intersection distance w s (i] ± , £ ) of the 
dual line through the triacontahedron. The index / denotes line densities. 

To obtain the average chain density for a given direction s one has to calculate the average 
intersection length w s . It is equivalent to the volume of the triacontahedron divided by the area 
^Tria,s of the triacontahedron projected onto the rj' 1 ('"'"-plane: 

VTria 

Pl, s OC W s = — . 



It is also possible to calculate frequencies of the chain densities. The pz jS -axis is divided up 
into intervals Spi yS . For each pi s the respective inverse image area in the ^(^-plane is calculated 
for the values Wg^ 1 - , f ) belonging to pi :S and 5pi yS . The total area is proportional to H{pi, s ). 
The chain intersection function w s {rj' L 1 Q^) is piecewise linear due to the flatness of the surface 
rhombs of the triacontahedron. It decreases monotonously in radial direction from the center of 
the interface (j}Tria,s to the boundaries. 

The frequency function may be derived analytically in the following way: The origins of the 
coordinate £ parallel to the symmetry axis s and of the coordinates r/^ and m the plane 
perpenticular to s are chosen in such a way that it coincides with the projection of the center 
of the triacontahedron. Then we compute the area of the interface F (^ ± ) of the triacontahedron 
with the same £"*- as a function of ^ for the half space ^ > 0. Due to the convexity of the 
triacontahedron the function -F(^" L ) decreases monotonously. From the inverse of the derivative 
of F^- 1 ) (Spi tS —>0) one obtains a function which is piecewise linear, has discontinuities and can 
be identified with the density function H(pi :S ) (Fig 



IV. ENERGY MEASURES 



A. Models of stability 



Several models currently exist to explain the stability of quasicrystals. In the deterministic 
energy model the internal energy U represents the thermodynamically stabilizing factor. Micro- 
scopic forces lead to matching rules ]l2| or overlapping cluster energies @,^9| that favour an ideal 
quasicrystalline tiling. In the non-deterministic entropic model stabilized by the entropy S no 
matching rules exist but the cells of the tiling do not leave any gaps and do not overlap. This 
is the random tiling model. In between is the semi-entropic model which is described by a free 
energy with contributions from internal energy and entropy: 

F(T) = U(T) - TS(T) . (2) 

The purely energetic model can be regarded as a low-temperature limit, the purely entropic model 
as a high-temperature limit of the semi-entropic model. 



B. Properties of energy measures 



The ideal quasilattice without any violation of matching rules represents the ground state of 
the energy model, taken at T=0. At finite temperatures thermally activated flips exist, mediating 
the transition between neighbouring states. The transition probability is given by the Boltzmann 
factor, which depends on the energy measure chosen. 

In this work we deal with the canonical random tiling ensembles. All configurations of an 
ensemble have the same volume and the same number of vertices due to a constant average 
slope of the de-Bruijn-hypersurface. The energy of the system is the ensemble average. The pure 
entropic random tiling model, on the other hand, is spezialized for microcanonical ensembles, since 
all configurations have the same energy. Starting from the occupation distribution of the energy 
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levels of a canonical random tiling ensemble in the thermodynamic equilibrium we can calculate 
the internal energy U(T) as the ensemble average U =< E > of the instantaneous energy E. 
The specific heat CV(T) can be derived from the variance of the occupied energy levels by the 
fluctuation-dissipation theorem 

8U 

< {5Ef >=< E 2 > - < E > 2 = k B T 2 C v = -jfi ■ 

The temperature-dependent internal energy and the specific heat are governed by the intrinsic 
properties of the system and the selected energy measure. The temperature dependence of the 
entropy density s(T) (precisely: entropy per quasilattice site[|) is given by the thermodynamical 
integration of cy : 

T 

s(T) = s v (T) + J dT' C -^f± + So . (3) 

o 

sy(T) is the entropy contribution of a change of volume. It vanishes for canonical random 
tilings, cy is the isochoric specific heat and sq:=s(T=0) the ground state entropy. At T— >oo we 
are in the limit of the pure random tiling model and we get the configuration entropy s(T=oo) 
=:Soo- The temperature variation of the specific heat depends on the energy measure, but not Soo- 
If one is interested only in the configuration entropy Soo, then the specific choice of the energy 
measure has no physical relevance if the first three of the following four conditions are fulfilled: 

1. The energy of any configuration is unique. 

2. The quasiperiodic reference tiling is a ground state. 

3. The energy measure is limited from above. 

4. The entropy of the ground state vanishes or is easy to calculate. 

The first condition requires an exact law to determine the energy of any microstate. The energy 
of a fixed microstate must not depend on how it is generated. Two configurations that differ only 
by a rigid translation have to be energetically identical. 

Especially global energy measures are endangered to violate the second condition and to render 
the quasiperiodic reference state unstable at T = 0. The third condition is a requirement for 

the integrability of CV jF^ as a function of T. If the energy measure is not limited, uncontrolled 
fluctuations of the energy in the high-temperature limit may exist. If they vary stronger than 
quadratic with T then 

< (SE) 2 > 
T 2 

diverges and with it cy as T—too. 

The last condition need not be met but is recommended for practical purposes. The ground 
state entropy so represents the number of states energetically equivalent to the quasiperiodic 
ground state. Since it is often complicated to calculate the ground state entropy it is useful to 
choose an energy measure where Sq vanishes or is easy to deduce. An example of a complicated 



1 Large letters indicate total quantities, small letters denote a quantity per vertex. For periodic approxi- 
mants with rhombohedral cells the number of cells and vertices is identical. 
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case of non-vanishing contribution so has been presented by Baake and Joseph p0[ for a two- 
dimensional octagonal quasicrystal. The locally defined energy measure is for the canonical 
vertices but positive for forbidden vertices and depends on the vertex type. Unfortunately there 
exist configurations which are not quasiperiodic, although they do not exhibit forbidden vertices. 
As a consequence 44% of the total entropy are contributed by the ground state. 

An example for a globally defined energy measure with an easily calculable ground state entropy 
is the ansatz of a harmonic oscillator potential in the orthogonal space as applied in this work. 
The energy is the sum of the squared distances of the dual quasilattice sites from th eir cen ter of 
mass. Up to a factor N it is equal to the variance in the orthogonal space (see Sec. Ill A ). The 
precise definition of the energy of the configuration a is 



E {a) = C 
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(4) 



The index denotes the ideal reference configuration, C=^p is an arbitrary normalization con- 
stant. 

This energy measure is called the harmonic energy measure. The variance for the ideal 
reference configuration is smaller than the variances for the overwhelming majority of the random 
tiling configurations a. But there is a tiny minority of configurations with a variance smaller 
than the value of the ideal tiling. Their atomic hypersurfaces are closer to a sphere than the 
triacontahedron. Such configurations may play a role at very low temperatures. To avoid energies 
less than the energy of the ideal tiling we have taken the absolute value in Eq. ||. For zero 
global phason strain in an ideal tiling this energy measure is not degenerate. But for periodic 
approximants the N possibilities (N: number of lattice points) to chose the origin of the unit cell 



yield a ground state entropy sq 



In TV 



per lattice point. It is obvious that sq vanishes in the 



thermodynamic limit. Strandburg j3[ has introduced a similar energy measure. But it was taken 
relative to a fixed point in the orthogonal space and not relative to the center of mass and therefore 
does not fulfil the criterium of finitencss of the energy measure (third condition). Without fixed 
boundaries of the system the whole distribution of the dual quasilattice points in the orthogonal 
space may drift and therefore yield a systematic contribution to the energy. 

Other global energy measures (not limited!) which we have used on a trial basis are the cubic 
energy measure 



and the quadratic energy measure 



N 



3=1 
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(the latter one has been introduced by Strandburg 0). M^a denotes the fixed center of mass 
of the triacontahedron in the orthogonal space. 

An example of a locally defined energy measure is the one that counts the violations of the 
alternation condition. This condition requires that along a worm two rhombohedra of the same 
type and orientation do not occur subsequently. This energy measure was used by Dotera and 
Steinhardt || and Gahler Q]. Other examples are the "simple energy model" which assigns the 
same energy to each forbidden verte, the "cluster energy model" which maximizes the frequency 
of certain favourite vertices and the "Tubinger mean-field model" . All three where used by Joseph 
in 0. 



7 



C. Pair interactions 



All energy measures discussed in Sec. [VB| contain a certain kind of arbitrariness since they do 



not depend on interatomic interactions but on the assignment of an increased energy value to ran- 
domized lattice configurations. Here we present test interactions which are motivated by pairwise 
interactions between atoms sitting on certain sites of the quasilattice (called "decoration" ) . We 
test whether the energetic equidistribution postulated in the random tiling model is justificable. 

The potential interactions used are the 12-6-Lennard- Jones potential modified by a cutoff func- 
tion to guaranty a smooth behaviour at the cutoff radius (r c = 4 in units of rhombohedra edges) . 
Another model is the potential of mean force 

v(r) = — ln(<?(r) + k), k > 0. 

g(r) is the radial density function and k an arbitrary parameter which screens the singularity 
caused by the zeros of g(r) (k = 100 in this work). To arrive at a potential which does not consist 
of purely delta-like minima due to the delta- like maxima of g(r) of an ideal quasicrystal one usually 
broadens the maxima with a gaussian. This yields an additional free parameter, the width (if all 
the maxima are broadened by the same function) 



The results of the test are presented in Sec. VIII C 1. It turns out that different random tiling 



samples are indeed equivalent with respect to the pair interactions, since the radial distribution 
function is nearly unchanged for the short distances (where the contibution to the potential energy 
would be strong), and at larger distances the pair potential is zero due to the cutoff distance. 



V. CONSTRUCTION OF RANDOM TILINGS 

A. Monte-Carlo method 

Simpleton flips directed in their frequency by the Monte-Carlo method allow the generation 
of equilibrium configurations at a given temperature. The randomization of the quasilattice is 
realized by a sequence of such flips which usually are associated with a change in energy. The flips 
play the role of elementary transition paths between different states of a system in contact with a 
heat bath. 

We are assuming that the method is ergodic [2f|. Although this has not yet been proven, there is 
no hint up to now of the opposite. The method has already been used extensively to generate and 
study equilibrium random tiling quasicrystals. It is the only one known that produces equilibrium 
ensembles. 

A simpleton flip does not generate or destroy tiles, it only rearranges them. Therefore the 
frequency of both types of rhombohedra in the tiling is remaining constant. The average orientation 
of the hypersurface and the global phason strain are not changed. 

During the Monte-Carlo simulation lattice points are selected at random. If a lattice point 
belonging to the simpleton vertex is hit, the energy of the original and the flipped configuration 
are compared. If the energy decreases, the flip is always carried out. If it increases, it is performed 
with a probability of exp(— ^MjO- 

During the flips the lattice points change their positions and the vertices their frequency. The 
simpleton vertex itself, for example, occurs with a frequency of 23.61% in the ideal tiling. This 
value decreases to a temperature dependent equilibrium value during the initial phase of a simu- 
lation. At T = oo the equilibrium frequency of the simpleton vertex is about 17.5%. 
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B. Boundary conditions 



Quasicrystals do not permit periodic boundary conditions. Hence one has to work with a finite 
patch of a quasicrystal. Then there are two possibilities to deal with the surface: either keep it 
fixed (fixed boundaries) or identify opposite sides (periodic approximants). 

In the case of fixed boundaries the surface lattice points are not allowed to move. In approx- 
imants all lattice points are mobile. Therefore the configuration entropy per lattice point with 
fixed boundaries is smaller than the entropy of approximants. The configuration entropy of the 
approximant, on the other hand, differs only slightly from the value of the ideal tiling. The fre- 
quency of vertices also is not much different between an ideal quasicrystal and an approximant 
tiling. The frequency of the simpleton, however, in patches with fixed boundaries may grow to 
up to 30% of that of the bulk vertices. For a discussion of the entropy dependence on boundary 
conditions see 22 . 

Periodic approximants show matching rule violations already in the ground state. They exhibit a 
periodic superstructure and for very small samples a remarkable deviation of the vertex frequencies. 
The intrinsic global phason strain changes the flip diffusion properties at low temperatures as the 
ground state phasons generate zero energy modes. The sheet magnetization (see Sec. [II C) in the 
case of a cubic approximant is useful only for the three two-fold directions parallel to the cubic 
cell axis. The acceptance region is no longer densely covered b ut has a lattice structure. This fact 



causes problems in the averaging procedures described in Sec. Ill E 



C. Addition method 



A new method to generate random tilings (especially non-equilibrium configurations) is the 
successive addition of tiles to the tube and slab structures within a finite region. 

For this method fixed boundaries are used. We start with a patch of a tiling with a fixed 
(possibly zero) global phason strain. This can be e.g. a piece of an ideal tiling or a single unit cell 
of a periodic approximant. All interior quasilattice points are removed. Then the volume is refilled 
in a systematic way without gaps and overlaps of the tiles. Since we start with an existing patch 
(microstate), it is certain that tilings of this surface must exist. The hollow volume represents the 
macrostate, the different possibilities to fill it are the microstates. 

The algorithm is started by selecting one five- fold direction (Fig. ||). First the ensemble of 
parallel slabs perpendicular to this direction has to be filled with rhombohedra. The borders of 
the slabs are given by a surface consisting of rhombi parallel to the five-fold axis chosen. A slab 
appears as a pentagonal Penrose rhombus pattern in projection. 

The surface of the hollow volume and the correlation between the slabs dictate the addition 
conditions. A violation sooner or later leads to gaps or overlaps, and the configuration has to be 
abandoned. 

We will now shortly sketch how the method works. The procedure can be broken up into two 
parts: filling the two-dimensional slabs, which is equivalent to performing the construction for 
generalized Penrose tilings, and then filling the space between the slabs. 

First of all we will explain how to proceed in two dimensions. In the Penrose tiling there are five 
sets of one-dimensional rows of rhombi with one type of edge in common. The average direction 
of the sets is perpendicular to the tiling vector directions. In the hollow volume we start with one 
of the directions and generate the rows one by one. The fixed surface determines the starting and 
terminating position of the rows (by the edges parallel to the selected tiling vector) . The difference 
vector (in higher dimensional coordinates) between initial and final point determines the number 
of tiles of each type in that row uniquely (by the coordinate components). We are only free to 
choose the sequential arrangement of the tiles. There are, however, restrictions: The second pair 
of edges of each rhombus must also fit into the configuration. This means that the difference 
vector between the rhombus and the respective edge on the surface of a neighbouring row must 
be fillable with rhombi. Otherwise, the configuration has to be abolished. After all the rows of 



9 



one direction have been created, the next direction is being dealt with until the whole patch is 
filled. The last direction is redundant, since the whole slab is filled as soon as the second-to-last 
direction is worked out. 

Up to a certain degree the process is running analogously in the three dimensional case. The 
slabs are being filled in the same way as for twodimensional Penrose-like rhombus tilings. We 
continue with the next slab perpendicular to the selected five-fold direction and generate slab by 
slab along this direction. Now we have to check the difference vectors of two pairs of edges to make 
sure that the space in between can be filled. But one must be careful as there are row directions 
crooked by oriented with respect to the row directions within each slab. For example, the row 
direction perpendicular to the rhombus spanned by the vectors t2 and t3 are crooked with respect 
to the corresponding direction of the vectors to and ti. To explain how to overcome this problem, 
however, takes more space than is available here. Interested readers are referred to [ p3[ . After all 
slabs of one direction are constructed the other directions are attacked until the whole volume is 
filled. Tilings generated by this method are called "addition tilings" . 

D. Comparison of the methods 

The acceptance domain for the ideal icosahedral rhombohedral tiling is a triacontahedron 
(Fig. ||) . In an equilibrium configuration generated by the simpleton flip method it is smeared out 
to an almost Gaussian distribution. Differences in the physical space are relatively small. There 
are some matching rule violations like pairs of adjacent rhombi of the same type, but the structure 
does not appear to be massively disordered. The phason fluctuations of the configurations are 
limited according to [|| . 

In a configuration generated by the addition method the deviation from quasipcriodicity is 
much stronger (Fig. |6|): There are microcrystallinc inclusions which indicate large- wavelength 
fluctuations in the higher-dimensional space. Many twin, domainlike or grain boundaries occur. 
The deviation from the ideal acceptance domain is strong and anisotropic. Configurations far from 
equilibrium are easy to construct. 

The efficiency of a specific method depends on its ergodicity properties. In a configuration 
generated by the addition method with microcrystalline inclusions there are few possibilities for 
flips. Very often the few flippablc lattice points only jump back and forth. There are bottle necks 
which prevent or strongly inhibit the structure from moving through phase space. 

Furthermore a configuration depends on its history. To generate independent configurations 
one has to wait until the autocorrelation is dropped off sufficiently. At high enough temperatures 
about 50 Monte-Carlo steps are sufficient. Fixed boundaries may increase the correlation times 
drastically. 

The anisotropic character of the configurations generated by the addition method is due to 
the fact that the slabs are filled sequentially. The deviation of the center of mass of the point 
distribution from the original acceptance domain depends on the sequence of construction of the 
slab and tube directions. On the other hand, it is only the sequential treatment that permits the 
tractability of the method. 

In the addition method the probabilities for adding a certain rhombohedron at a certain place 
have not been taken into account. Therefore, the configurations generated by this method in a 
selected hollow volume can not be considered as the sample survey of a equilibrium ensemble. To 
calculate the probabilities would be impossible since the number of realizations is prohibitively 
high. 

The addition method allows the generation of configurations that are never reached in acceptable 
times by the simpleton flip process. The method could be used to generate random tilings with 
constraints. 

Experimentally a number of decagonal quasicrystals have been found that maybe could be 
described by configurations generated by the addition method p^-|2^]. The distribution of the 
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quasilattice points in the orthogonal space is not isotropic as in the samples produced by the 
simpleton method but strongly anisotropic. 



VI. RANDOM TILING TRANSITION AND STRUCTURE FUNCTIONS 



Signals for a random tiling transition can be obtained in a number of ways. The thermodynamic 
functions u, cy and s may show a characteristic behaviour. For example a plot of the specific heat 
cy versus temperature for different sample sizes (finite-size scaling) may indicate the possibility 
of a random tiling transition. The specific heat cy should diverge in case of a second order phase 
transition. 

For any desired temperature a histogram of frequency of the energy values can be calculated 
during the simulation. The histogram for infinite temperature yields an approximation for the 
density of energy states since at infinite temperatures all energy levels should be occupied with 
equal probability. Therefore it may be used to explain the behaviour of the specific heat as a 
function of T (see also Sec. [VIII Bl[ ). 

In addition to the specific characterization criteria for quasicrystals described in Sec. HI there 
are more general functions that can indicate a phase transition. We have applied the following: 



1. Self diffusion coefficient 



The average squared difference ((r(i) — r(0)) 2 ) of the coordinates of all quasilattice points from 
their initial position is expected to depend linearly on time at long enough simulation times. The 
slope of the function is proportional to the self-diffusion coefficient D. A change of the slope may 
indicate a phase transition which involves the change of energy barriers for example. 



2. The Binder order parameter 

The Binder order parameter [§§ SW(T) for a given sample size N and temperature T is 
defined by 

bW(t):=1 _ Mf\T) 



3(M 2 W(T)) 

where M 2 (A °(T) and m[ N \t) are the second and fourth moment 

M[ N \T)= /d M ( P ( M )W(T)y 



of the probability p(n) of a microstate of the sheet magnetization ji belonging to a macrostate 
with sheet magnetization M. The Binder order parameter B^ N >(T) is plotted as a function of T 
and yields a set of curves parametrized by N. A unique intersection point of the curves points 
to a second order phase transition and yields the transition temperature. If there is no clear 
intersection point only an interval for a possible phase transition can be determined. 



A. Structure functions 
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1. Radial distribution function G(r) 

The radial distribution function G(r) of lattice points or atoms (in case of a decorated quasi- 
lattice) consists of discrete delta- type maxima and is related to the radial density function g{r) 
by G(r) = 4irr 2 g(r). The Fourier transform of the radially averaged G(r) is the radial structure 
factor I(k). The potential energy u of an quasicrystal can be directly calculated from G(r) if there 
are only pair interactions v(r): 

R 

u = Y J G{r)v{r) (5) 

r=0 

where R is the cutoff radius of the potential. If there are different types of atoms, then additional 
sums over the different types of interaction pairs must be taken. 

2. Structure factor intensity I(k) 

Flips of the tiles modify the occupation density of atomic chains and atomic layers. Diffraction 
diagrams are dominated by the long-range translational and orientational ordering of quasicrystals. 
Therefore an alterations of the quasiperiodicity changes the width and intensity of reflexions. 

Phason fluctuations furthermore generate diffuse scattering. Due to the limitation of the phason 
fluctuations Q the maximal intensities are barely influenced. Diffuse scattering is governed by 
the phason-elastic constants. 

3. Diffraction diagrams 

Diffraction diagrams for random tilings were calculated for randomly selected sample configura- 
tions filed during the simulation process at a given temperature T and size N and averaged over 
the whole set chosen. The fluctuations between the selected configurations turned out to be small 
at a given T and N. 

VII. SIMULATION PROCESS AND MODELS STUDIED 

The internal energy u, the specific heat cy, the entropy s as well as the sheet magnetization 
M, the corresponding susceptibility \ an d the self-diffusion coefficient D, all dependent on the 
temperature T, are thermodynamical state variables of a random tiling ensemble. The structure 
of a single tiling, but also the whole ensemble can be characterized by structural functions. In this 
work we have used the diffraction-intensity I(k) (without paying attention to diffuse line shapes), 
the relative frequencies of vertex environments, the axial and planar structures, the variance Q 
(phason fluctuations) and the radial distribution function G(r). The quantities u, M, Q and 
{(r(t) — r(0)) 2 ) ^] are sampled by registring their values at regular intervals of length At during 
the simulation. From the collected data all above mentioned state variables are calculated. Data 
of microconfigurations are stored every 1000-Ai Monte-Carlo step. The results may deviate from 
the thermodynamical limit due to finite-size-effects. If the sample size grows, the values should 
converge to the thermodynamic limit. 

An initial configuration (e.g. an ideal Ammann-Kramer-Penrose approximant or an addition 
configuration) was equilibrated over a typical thermalization time of 10 3 to 10 steps before the real 



2 The simulation time is set to zero after thermalization. 
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simulation was started. The thermalization time was chosen according to the saturation behavior 
of the energy, for example. At low temperatures longer equilibration and simulation times and 
sampling intervals were used. The length At was checked by the decay of the autocorrelation 
functions. For high temperatures At was set typically about 50 Monte-Carlo moves per lattice 
point, while for lower temperatures 200 moves were taken due to temporal correlations. The entire 
duration of a simulation run at given temperature typically has been about 25000-Ai to 30000-Ai 
simulation steps. 



A. System sizes and interaction types 



The following systems have been studied in our simulations: 

1. A finite patch of a quasiperiodic lattice with 4403 lattice points (3507 cells), fixed boundaries 
and the cubic energy measure. Starting configuration has been the ideal quasilattice. 

2. A cubic approximant with 2440 lattice points, periodic boundaries and the quadratic energy 
measure. Starting configuration was the ideal approximant. 

3. Two five-fold approximants with 890 and 1440 lattice points, periodic boundaries and the 
harmonic energy measure. Starting configuration also has been the ideal approximant. 

4. Five cubic approximants with 136/712 (n=3), 576/3016 (n=4), 2440/12776 (n=5), 
10336/54120 (n=6) and 43784/229256 (n=7) lattice points/atoms in the binary decoration 
modeF 1 29 1 (n is the generation), periodic boundaries and the harmonic energy measure. In 



this category we used ideal approximants as well as addition tilings as starting configura- 
tions. Here we also calculated the self-diffusion coefficient D(T), the magnetization M(T) 
and the susceptibility x(T)- The results are independent of the starting condition in the 
case of periodic boundaries: Ideal quasicrystals and addition tilings yield the same ensemble 
averages^, proving that the system had been well equilibrated before data collection. 



VIII. RESULTS 



A. Direct calculation of the configurational entropy 

For small clusters of rhombohedra with fixed surface it is possible to calculate the configurational 
entropy s directly, which we have performed for the five smallest polyhcdra with triacontahedral 
shape and edge lengths of one or two. The results are presented in Tab. |[ 

The clusters are too small to permit a reasonable extrapolation to the thermodynamic limit. 
The value of the entropy is also too small due to the fixed boundaries. 



B. Thermodynamic functions 



3 In the binary model each vertex and edge center is decorated with a small atom, and the long diagonal 
of the prolate rhombohedron is decorated with two large atoms 
4 which is not true in the case of fixed boundaries 
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1. Internal energy, specific heat, entropy 



In sample one the ground state entropy vanishes. The internal energy grows with T and saturates 
at high T. Due the fixed boundaries the saturation limit is finite although the energy measure 
is unlimited. The specific heat grows until T w 1. Then it decreases and goes to zero at high T 
i. e. we observe a Schottky anomaly as is typical for two level systems. The entropy at infinite 
temperature approaches Soo » 0.1194 ± 0.015. 

Sample two behaves similar, but the saturation of the internal energy for T > 10 is not so clearly 
visible and may indeed not occur since the energy measure is not limited and periodic boundaries 
were applied in this case. The decay of the specific heat is slower which is also a consequence of 
the energy measure. The entropy at infinite temperature approaches m 0.2621 ± 0.015. 

Samples three and four are again quite similar. The internal energy grows monotonously, and 
a clear saturation becomes visible (Fig. fj]) . The saturation value of u depends on the size of the 
sample. The limit value can be derived from limits of the variance: Q(T=oo) = 1.73 ± 0.01 || 
leads to u between 1.97 and 2.05. A Schottky anomaly in the cy-Plot again is present but the 
specific heat shows an additional bump above the maximum (Fig.g). Such a behaviour is known 
for few-level systems (for example three-level systems) with sufficiently separated levels. We have 
mapped the distribution of the energy levels for T—oo. No indication of discrete energy levels 
was found, only an asymmetry of the distribution with a smaller slope at higher energies could 
be observed. There are other energy measures which exhibit no visible asymmetry in the 
energy distribution and no bump in the specific heat. This, however, is in our opinion not a clear 
explanation for the phenomenon. The value of the maximum of cy is not significantly dependent 
on sample size. The increase in Sqo is caused only by the growing width of maximum (Fig. |^). 
Entropy values at T = oo are listed in Tab. [|. Since no divergence occurs with increasing sample 
size, we have no hint for a second order phase transition. 

2. Simpleton magnetization, simpleton susceptibility, Binder order parameter 

The temperature dependence of the energy fluctuations and the specific heat do not give a clear 
indication for a phase transition. There is no divergence of the maximum of cy. The reason may 
be that the intrinsic divergence of the specific heat with sample size, if any, is very weak. For 
further insight we calculated the sheet magnetization M and the susceptibility x since the latter 
shows a much more pronounced divergence behaviour. 

In the samples of class four the magnetization saturates at high T, and the minimum value 
decreases with size (Fig. |l(|). At higher temperatures, between T = 3 and 10, it increases again. 
This behavior we attribute to the finite size of the samples. The value of the maximum of the 
susceptibility grows almost linearly with the generation n and moves to lower T (Fig. [ll|). It is 
not yet clear if the relation Xmax{n) °c (n — tiq) is valid for n > 7. If yes, this would be a slow 

divergence (more precisely: Xmax(N) is about proportional to yN where N is the number of 

lattice cells and r— ^ > i^~ - ) ■ With our current simulation programs it is not possible to deal with 
n > 8. 

The Binder order parameter does not lead to a decision either, since there is no clear intersection 
point visible (Fig. |l2|). The trend of B^ N '(T) changes at n = 5, similar to the magnetization 
behaviour. The most probable intersection point is between T = 0.3 and 0.4 if it exists. 

But the tendency of Xmax with n and the behaviour of B^ N \T) indicate that a phase transition 
is much more likely to occur at T — 0. 
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3. Self-diffusion coefficient 



The mean square displacement (ir{t) — r(0)) 2 ) grows linearly with time t, indicating a normal 
diffusion behaviour and allowing the calculation of the diffusion coefficient D. At temperatures 
lower than T w 1 the extrapolation is difficult due to large fluctuations. D becomes unmeasurably 
small below T = 0.5. The diffusion coefficient forms a plateau at T ^ 1 for n = 4, 5, 6 in the 
Arrhenius plot (Fig. [[3]). There may be several reasons for this behaviour: 

First, there are energy barriers between different tiling configurations due to the harmonic 
energy measure, which at low temperatures lower the mobility of lattice points for flips of higher 
energy. In the range of the plateau the probability for a flip only occasionally suffices to overcome 
the barriers which play no role at high temperatures. 

Second, a phase transition may occur which changes the slope in the Arrhenius plot. Gahlcr 
has also observed a change of the slope in case of the energy measure with the alternation condition 
which turns out to be a phase transition since other response functions like the susceptibility and 
the specific heat definitely yielded a divergence at the same temperature. 

Last, there is an explanation which comes from the fact that the number of flippable lattice 
points (number of simpletons) changes with temperature. We can distinguish four ranges: 

1. The number of simpletons is about 23% in the range < T J$ 1. 

2. In the range 1 J$ T ;$ 10 we find a nearly logarithmic decrease of the number of simpletons. 
The plateau of D(T) is clearly seen here. 

3. In the range 10 J$ T ;$ 100 the approach of the frequency of the simpleton to a constant 
value leads to an increase of the negative slope of D(T). 

4. Above T ^ 100 the number of simpletons is constant at w 17.5%. 

In the approximants the behaviour of D(T) is obscured to some degree by zero energy modes 
caused by periodic boundaries. These modes become less and less important at larger sizes, but 
suppress the plateau for small sample sizes. 

Within the framework of the random tiling model we have studied only flip diffusion. No other 
diffusion mechanism, in particular vacancy diffusion, can be introduced in this way. The latter 
mechanism is expected to be the dominant diffusion process at least above 600°C, but there are 
indications that the new mechanism plays a role below that temperature [ ^0[ . 

C. Structure functions 

1. Radial density function, radial structure factor, pair interaction energy 

Structure functions have been calculated for samples of type four and size n = 6, containing 
10336 lattice points, for the primitive monoatomic and the binary decoration [ p9| . There are only 
small changes in the radial density function g(r). The largest differences of the order of 12% occur 
for bonds between two large atoms, but they are the least significant ones for stability due to their 
small number. The energetically most important nearest-neighbour bonds change even less in all 
cases, since a large portion of them lives within the rhombohedra and thus is only transferred to 
another position by a flip, but the frequency itself does not change. The differences we observe 
in the addition tiling are somewhat larger, but still very small. Figure |lj shows the results for a 
monatomic sample of size n = 5. 

The trend of the radial structure factor I(k) which corresponds to the FouRiER-transformation 
of g(r) is similar to that of the radial density function. Due the smallness of the sample there are 
rather large finite size effects which do not permit a quantitative comparison. 

We further observe a very small dependence of the pair energy function (see Eq. |^) on size and 
T due to the relation between the pair energy and the pair distribution function. It is interesting 
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to note that the potential energy of the random tiling ensemble is at T = oo about 0.4% lower 
than for the ideal quasicrystal for a large class of Lennard- JONES-like potentials, whereas the 
fluctuations within the random tiling ensemble are of the order of 0.1%. Therefore simple pair 
potentials obviously favour energetical equidistribution, i. e. the random tiling model. 



2. Variance, vertex statistics 



The variance behaves roughly linear with respect to u which is obvious if one compares their 
definition in Eqs. ^ and ^. The variance therefore shows a similar behaviour as the internal energy. 

The change in the statistics of the vertices with temperature is much more important, as we 
have s een in th e discussion of the diffusion behaviour and in the role of the simpleton frequency 
in Sec VIII B 3 . The changes of the vertex frequencies (Notation as in |3l| ) can be summed up as 



follows: 

1. The simpleton vertex 1 (452)^ decreases from 23.61% to 17.5%. 

2. Vertex 3 (670) decreases from 23.61 % to less than 10%. 

3. Vertex 2 (561) increases from 23.61 % to 29 % at T = 8 and then decreases to 26%. 

4. The sum of frequencies of vertex 1 and vertex 2 is roughly constant. 

5. The forbidden vertices (all together) increase from 0% to more then 40%. 

6. The vertex with the highest symmetry where 20 prolate rhombohedra meet (vertex 24, in |3l]| 
also called twelve- fold sites) decreases in frequency from 1.2% close to extinction. 

The addition tilings show a strong deviation from the ideal distribution. For the sizes n = 4, 5, 6 
we obtain for the simpleton frequencies: 16.7 %, 15.0 %, 12.2 %, and for the frequencies of the 
forbidden vertices: 36.3 %, 43.2 % and 55.3%, respectively. 

We note that the statistical error of the vertex statistic is — in the case of equilibrium random 
tilings — largely independent of the sample size and already for n>4 very small. 



3. Plane and axial structures 



Plane and axial structures have been calculated for a sample of type four with the size n = 6, 
containing 10336 lattice points. The area function A s has been calculated from the dual 
lattice points in the orthogonal space, and the hist ogram of the densities H(pp) has been derived 
by numerical differentiation as described in Sec. [II D| . The analysis has been carried out at 



temperatures T = 0.3, 1,3,... 1000 for planes perpendicular to two-, three- and fivefold planes. 

There is a continuous transition from an area function A s (£ ) which is dependent on the sym- 
metry directions to an almost isotropic gaussian shape at high temperature, indicating that the 
acceptance domain is smeared out isotropically. At T = oo the histogram of densities H{pp) also 
becomes almost independent of the symmetry direction. H(pp) is strongly peaked at very low 
and very high plane densities which indicates that very weakly and very densely occupied planes 
dominate. 



5 The first digit is the number of vertices connected by an edge, the second the number of vertices 
connected by a short face diagonal and the last the number of atoms connected by a short diagonal of the 
oblate rhombohedron 
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In the five-fold direction we observe the smallest changes. Strong differences between the analyt- 
ical and numerical functions, especially at poles and singularities, for three- and two-fold directions 
prohibit quantitative predictions. 

The addition tilings show a broader distribution in A s (£ ) and an anisotropy. Additional minor 
maxima occur in H(pp), but they alternate from one realization to the other. The distribution of 
dual lattice points of one tiling is anisotropic by construction. To obtain A s (£-*-) we have averaged 
over all symmetry-equivalent directions. 

We could not find any useful results for axial structures. The size of the samples appears to be 
too small to give a satisfying resolution for histograms. Strong differences between analytical and 
numerical H(pi) functions for T=0 also indicate that the results are not reliable. 

4- Diffraction patterns 

Diffraction patterns (Bragg scattering without observation of diffuse line shapes which are char- 
acteristic for random tilings) have been calculated for a sample of type four with the size n = 5, 
containing 2440 lattice points. The analysis has been carried out at temperatures T — 0.3, 1, 3, 
10, and at T = 10000 for planes perpendicular to two-, three- and fivefold planes. The intensities 
are only weakly dependent on temperature as the acceptance domain is only weakly smeared out 
even at high temperatures. This result is consistent with Tangs observation of limited phason 
fluctuations in three dimensions 

In the case of the addition tilings stronger differences exist: Some of the reflexions are weaker, 
and new ones arise. This indicates that there is already a large phasonic disorder in the addition 
tilings which generates an extended cloud of dual lattice points. 

IX. CONCLUSIONS 

We have investigated a geometric and thermodynamic model of disordered quasilattices, which 
could explain the structure of icosahedral metal alloys. 

In this random tiling model there are no matching rules, in contrast to the ordered quasiperiodic 
structures. The only requirement is that the physical space must be filled with rigid tiles without 
gaps or overlaps. 

The results are summarized as follows: 

1. The internal energy u increases and saturates at high temperatures. Moo converges towards 
a limit for N— >oo, which can be estimated from values of the phason fluctuation variance 
for the pure entropic random tiling ensemble. The entropy obtained from the specific 
heat cy by thermodynamic integration is about 0.24±0.01, in agreement with results from 
many other sources |3p|Jl]|. The specific heat cy shows a Schottky anomaly. It is not 
clear whether it diverges. The increase of the entropy with the system size comes from the 
broadening of the maximum. 

2. The sheet magnetization decreases to a minimum with the temperature and slightly grows 
to a saturation value. In the thermodynamic limit it should vanish without an intermediate 
minimum. The susceptibility diverges slowly, the maximum shifts to smaller temperatures. 
The Binder order parameter does not exhibit a unique intersection point. We assume that 
the transition temperature is finite, hence the random tiling transition could be analogous 
to the transition in a two-dimensional Ising ferromagnet. 

3. The self-diffusion coefficient displays a plateau in the central temperature range indicating 
energy barriers for certain flips due to the harmonic energy measure. On the other hand there 
exist correlations between the temperature dependence of the self-diffusion coefficient and 
the frequencies of simpletons per lattice point. The Arrhenius plot is deviating strongly from 
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that assumed by Kalugin and Katz pj] : Where these authors are plotting a steep increase, we 
are observing the plateau. There are zero energy modes due to periodic boundary conditions, 
which lead to a suppression of the plateau for small approximants. 

4. Henleys postulate of the finiteness of the phason fluctuations is demonstrated by the Monte- 
Carlo-simpleton-flip tilings. The addition tilings show larger variances and deviations of the 
vertex frequencies from the equilibrium values. There are also changes in the diffraction 
patterns caused by phason fluctuations. 

5. The frequencies of the densities of planes in the equilibrium tilings are temperature dependent 
and very weakly and very densely decorated planes dominate. Addition tilings show further 
minor maxima, dependent on the realization. 

6. Radial structure functions depend only weakly on the configuration. This is due to the 
rigidness of the cells. Pair interactions realize an equidistribution of all configurations — a 
possible realization of the pure entropic random tiling model. 

The behaviour of cy and \ obviously depends strongly on the type of energy measure used. The 
alternation condition seems to work much better |ll| — maybe as consequence of its closer simi- 
larity with Ising interaction models in comparison with the harmonic energy measure. However, 
for the simple energy model it was also not possible to decide if a phase transition occurs, 
since no divergence of the specific heat could be observed. 
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TABLE I. Configurational entropy for triacontahedral clusters with fixed surface. The first column 
codes the type of polyhedron by the length of edges in units of the six icosahedral basis vectors. The 
second column gives the number of configurations, the third the number of rhombohedra. The last 
column contains the configurational entropy per lattice point. 



edge lengths 


# configs. 


# rhombo. 


s 


111111 


160 


20 


0.2538 


211111 


1280 


30 


0.2385 


112211 


22381 


44 


0.2276 


121121 


22981 


44 


0.2282 


221121 


1268131 


63 


0.2231 



TABLE II. Some thermodynamic parameters for the harmonic interaction model with periodic boun- 
dary conditions. The first column provides the generation for cubic approximants, the second column the 
number of lattice points, the following columns list the entropy s x and the variance f2(T=oo). The rows 
"890" and "1440" denote pentagonal approximants. 



n 


size 


Soo 


n(T = oo) 


3 


136 


0.1816878 


1.475 ± 0.015 


4 


576 


0.2032964 


1.56 ± 0.01 




890 


0.2086351 


1.60± 0.02 




1440 


0.2076728 


1.615± 0.02 


5 


2440 


0.2114427 


1.635 ± 0.01 


6 


10336 


0.2349887 


1.675 ± 0.01 


7 


43784 


0.2367582 


1.70 ± 0.01 




limit 


0.24±0.01 


1.73 ± 0.01 
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FIG. 1. Simpleton flip 



FIG. 2. Worm of rhombohedra 
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FIG. 3. Plane structures in the ideal tiling. The left column contains the functions v4 s (£ x ), the right 
column displays the histograms of plane densities H{pF, s )- 
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FIG. 4. Chain structures in the ideal tiling. The histograms of chain densities H(pi tS ) are given. 
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FIG. 5. Intermediate state in the addition process. The star indicates the tiling vector directions. 
The double arrow indicates the direction which has already been filled. The small arrows indicate the 
correlations to the direction which will be filled next. 
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FIG. 6. Comparison between different tilings: a) perfect quasicrystal, b) generated by flip process, c) 
generated by addition process. The corresponding acceptance domains are labeled by d), e) and f). 
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FIG. 7. Internal energy. Sizes are indicated by the number of lattice points. 



FIG. 8. Specific heat. Sizes are indicated by the number of lattice points. 
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FIG. 9. Configurational entropy. Sizes are indicated by the number of lattice points. 



FIG. 10. Magnetization. Sizes are indicated by the number of lattice points. 
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FIG. 11. Susceptibility. Sizes are indicated by the number of lattice points. 



FIG. 12. Binder order parameter. Sizes are indicated by the number of lattice points and the generation 
parameters n. 
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FIG. 13. Self-diffusion coefficient. Sizes are indicated by the number of lattice points. 



FfG. 14. Radial density function (RDF) of the rhombohedra lattice points. The peaks represent the 
values for the ideal tiling, the crosses denote the difference for the random tiling at T = oo. 
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